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Tunnelling between non-centrosymmetric superconductors with significant spin-orbit splitting is 
studied theoretically in a two-band treatment of the problem. We find that the critical Josephson 
current may be modulated by changing the relative angle between the vectors describing absence of 
inversion symmetry on each side of the junction. The presence of two gaps also results in multiple 
steps in the quasiparticle current-voltage characteristics. We argue that both these effects may 
help to determine the pairing states in materials like CePtaSi, Ulr and Cd2Re207. We propose 
experimental tests of these ideas, including scanning tunnelling microscopy. 
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Superconductors where inversion symmetry is absent 
and spin-orbit splitting is significant have recently at- 
tracted considerable attentioniSiiMiiiiSiLli Much of this 
interest was initiated by the discovery of superconduc- 
tivity in CePtaSi^ and UlxpiS In addition to not hav- 
ing an inversion center, band structure calculations on 
CePtaSi* have shown that spin-orbit coupling splits oth- 
erwise degenerate bands by 50-200 meV near the Fermi 
level. This is much larger than k^T c , where T c is the 
critical temperature for superconductivity. This must be 
taken into account when describing superconductivity in 
these systems. These materials also order magnetically, 
which could influence the nature of the superconduct- 
ing state. However, at least for CcPtaSi, there seems to 
be little communication between superconductivity and 
magnetic order >ii 

Another superconductor of interest in this context is 
Cd2Re207. This material has a structural phase transi- 
tion where a center of symmetry is lost^* When certain 
ions in the unit cell are displaced throughout the lattice, 
internal electric fields are induced, giving rise to spin- 
orbit splitting of spin-degenerate states. Calculations 
and photocmission studies^* have indicated that this 
splitting will have a significant influence on the electronic 
band structure. Thus, Cd2Re207 is similar to the mate- 
rials mentioned above, although simpler, since it shows 
no sign of magnetic order^a Another pyrochlore super- 
conductor that might fall into this category is KOS2O6. 15 

The spin-orbit splitting of otherwise degenerate bands 
demands a two-band description of superconductiv- 
ity in these materials. An exotic feature of non- 
centrosymmetric superconductors with large spin-orbit 
splitting is the possible absence of a definite parity of the 
superconducting statei*£ Experiments have indicated 
that CePtsSi might be in such a pairing state, a linear 
combination of spin-singlet and spin-triplet states, and 
that the gap may contain line nodes£ Cd2Re207 and 
KOS2O6 seems to be nodeless, however i^ilk 

In the present study, we will investigate tunnelling cur- 
rents between two superconductors where a two-band de- 
scription is necessary in both systems. Junctions involv- 



ing one superconductor with spin-orbit split bands have 
been studied in Refs. and^| We will restrict ourselves 
to intraband Cooper pairing without specifying the mi- 
croscopic mechanism responsible for this. See Refs. ^3 
for related work on MgP>2 junctions. Our main focus will 
be on non-centrosymmetric superconductors with spin- 
orbit split bands and we will specialise to this case when 
needed. We find that the critical Josephson current may 
be modulated by changing the angle between the vectors 
describing absence of inversion symmetry on each side. 
This effect is analogous to tunnelling magnetoresistance 
in ferromagnetic tunnel j unctions We also calculate the 
quasiparticle current. For temperatures close to zero, the 
current- voltage diagram may contain several discontinu- 
ities determined by the relative size of the two gaps. We 
claim that both these results may help to determine the 
possibly novel properties of the superconducting state in 
materials like CePtaSi, Ulr and Cd2Re207. 

The Hamiltonian considered is H = + Hgc, where 
i?sc describes superconductivity. The normal state 
Hamiltonian is 

H N = ^2<pl[e k + B k -(T}4> k , (1) 

k 

where (f> k = (c k p c^), e k is the band dispersion, and 
the vector er consists of the three Pauli matrices. We 
name the spin quantisation axis the z-axis. The vector 
Bfc removes the spin degeneracy from the band e k . By 
a transformation to a basis 4> k = (ct k , c_ k ) where 
is diagonal, one finds H N = J2\=±, k £\,kc\ k c\ yk . The 
quasiparticle spectrum is E± tk — e k ± \B k \. We define 
B k ,± = B ktX ± iB k , y = B k ^±e ±lipk . 

The vector B k has the property B k — —B kl - where 
B k characterises the absence of inversion symmetry in 
the crystal. The origin may be that ions are removed 
from high-symmetry positions, as in ferroelectrics^* 
leading to internal electric fields and thus increased spin- 
orbit coupling. 19 20 To establish the form of B k , point 
group symmetry considerations may be employed* and 
B k will depend on the direction in which the ions are 
displaced. 
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An electron with momentum k will align its spin par- 
allel or antiparallel to B k . In a free electron model with 
the Rashba interaction^ the ID density of states for the 
+ and — band at the Fermi level are equaliSi Still, we 
allow these to be unequal, which is the general case. 

Let us now turn to the term responsible for supercon- 
ductivity, -ffsc- We write down the interaction in terms 
of the long-lived excitations in the normal state 



Xfi,kk' 



This model contains only intraband Cooper pairing. In- 
terband Cooper pairs are strongly suppressed if the 
spin-orbit splitting is much larger than the supercon- 
ducting gaps, even though the two bands may touch 
at some isolated points on the Fermi surfaced This 
is the limit we are investigating. Defining Aa,& = 
~X^,fc' v ^( k y fc'Xc/^k'C^-fc'), the standard mean field 
approach gives the total Hamiltonian 



H = 



\,k 



£\,kc\ k cx,k + 



A,fc°A,-fc 



Note that A Aj _fc = — A A .fc follows from the fermionic an- 
ticommutation relations. In Equation ©, the two bands 
are decoupled, resulting in Green's functions diagonal in 
the band indices. This is a result of the mean field ap- 
proximation. A±_ k are in general not independent, but 
related through the self-consistency equations due to the 
possibility of interband pair scattering, - 

The relation = — B k ensures that states with op- 
posite momenta within a band have opposite spins. For a 
spin-1/2 state, the time reversal operator is JC = —io y K,Q, 
where /Co denotes complex conjugation. Let the original 
operators transform according to K, : c\, a = — crc^_ k _ a 
under time reversal. The effect of time reversal on the 
new operators then becomes K. : ct h = i A ,fcC A _ k , where 

t\,k = e~ XiVh . This means that if \\,k is the order 
parameter for pairs of time reversed states, one finds 
A A ,fc = t x ,kX\,k- This gives \\,k = X\-k- Thus, \\,k 
may be expanded in terms of even basis functions of ir- 
reducible representations of the space group£ 

Define the matrix A k whose elements are the gap func- 
tions Afc !(Tcr / in the original basis where spin is quantised 
along the z-axis. This may be written as 



A fc = ?7fc. s g + rj k ^ T (B k ■ er) g, 



(4) 



where g = —\o y . The first term is symmetric in mo- 
mentum space and antisymmetric in spin space, whereas 
the opposite is the case for the last term. Thus, in 
the absence of spatial inversion symmetry, the order pa- 
rameters in a spin basis have no definite parity, but 
is in general a linear combination of singlet (S) and 
triplet (T) partsi2i& The singlet and triplet compo- 
nents are determined by r]k,s = (x+,fe + X-,fe) /2 and 
Vk,T = (x+.k — X-,k) /2, respectively. This means that 



knowledge of x±.k and B k could help determine the gap 
structure and the symmetry of the pairing state. For non- 
centrosymmetric materials like CePtaSi, this is currently 
a matter of intense studyA^SiS 

The normal and anomalous Green's functions for each 
band are G x (k,uj n ) = -(iw„ + 6v,fc)/(^ + £i\fc + IXA,fc| 2 ) 
and T\(k,uj n ) = t x ,kX\,k/{^l + £ A)fe + IXA,fe| 2 ), respec- 
tively. These are defined in the standard way, see e.g. 
Ref. H uj n is a fermion Matsubara frequency, ^\,k = 
£\,k — M an d M is the chemical potential. 

Consider tunnelling between two such superconduc- 
tors, A and B. Let system A be described by the Hamilto- 
nian ©. The Hamiltonian of system B is defined equiv- 
alently, only with Cka,c\ k — ► d pai d PtP . Moreover, we 
allow B k and B p to be different. Consequently, even 
if fe = p, the spin in a state + or — may be different 
on sides A and B. The tunnelling Hamiltonian is ifr = 

Yl<\pMp {^kp^XM^P.P 
Ap 



J fcp u p : p L A,fc 



The tunnelling 



matrix elements T kp depends strongly on the direction 
of k and p. 22 Tunnelling is much more probable for a 
momentum normal to the interface rather than parallel 
to iti22iSi24 If we assume that spin is conserved in the 
(3) tunnelling process, i.e. H T = Yl,kp,a T k v ^ ka A V a + h.c, 



we find IT, 



Api2 



kp 



\T kp \ 2 (l + \pB£-B*)/2. 



The number operator for band A in system A is N^- = 
Efccl.fcCA.fe- We define TVj = i[H T ,N^], such that the 

charge current is I(t) = —eJ2x(^D- To lowest order 
in the tunnelling matrix elements, standard theory gives 
I(t) = J qp + Ij(t), where I qp = -2e £ A Im and 
Ij(t) = 2e£ A Im [^\(eV)e 2ieVt ]. The voltage is eV = 
Ma — Mb ■ In the Matsubara formalism, we have 

**K) = \ £ |f fc A ;i 2 ^(fc,^-^)c? p B (p,c n ), 

P k P (5) 

p,U) n 



kp 



where uj v — > eV + i + is a boson Matsubara frequency. 
We have assumed that the bulk Green's functions may 
be used, neglecting boundary effects. Such effects could 
however be of importance in these systems^ due to the 
possibility of subgap surface bound states or distortions 
of the order parameters close to the surface. 

Time-reversal symmetry of ifp gives T_ k p = 

Tk£* t\,ktp?p^ These phase factors will cancel the ones 
from the anomalous Green's functions in Equation JSJl, 
which shows that each band A may behave as a singlet 
superconductor with gap function x\,k~^ 

We take the continuum limit^i and assume that 
A/" A (£, f2), the angle-resolved density of states in band A 
in the non-superconducting phase, is constant. 

The gap xa(£, ^) depends on both energy and the di- 
rection of momentum. Neglecting the energy dependence 



3 



is standardised The tunnelling matrix elements ensures 
that momenta approximately perpendicular to the inter- 
face will dominate d 2 ' 2 ?' 24 We therefore let Xa(£j ^) ~ Xa, 
the value at the Fermi level and directions normal to the 
interface (remember that XA.fe = XA,-fc)- This is exact 
if the gaps are isotropic or if the tunnelling is strictly 
one-dimensional. It could also be a good approximation 
if the variations of x+ and X- are small in the region 
around normal incidence. We define xa = Ixa |e 11?A . 

The energy dependence of the tunnelling matrix ele- 
ments may be neglected. We will need the quantity t\ p = 
J dQ A J dfl B \f x P(fl A , fl B )\ 2 . Let us look at a specific ex- 
ample, where |Tfe p | 2 ~ \T\ 2 k±p±Q(k±p±)m— In addition, 
we choose the Rashba interaction B^ = a(n A x fe)>iS We 
let h A and n B , and consequently the nodes of B^ and 
.Bp, point parallel to the interface. The Rashba inter- 
action appears to be an appropriate choice for CePtaS: 4 
and Cd2Rc2C>7. 5 Define the angle ( by cos£ = n A ■ n B . 
This gives 



\T\ 

t\ p = -y- (i + x V COS c) , 



(6) 



with x w 0.6d Numerical integration indicates that 
Equation JBJ is a very good approximation also when 
parallel momentum is conserved^ In general, it seems 
reasonable that if mostly electrons near normal incidence 
contribute to the current, Equation © is applicable with 
cos C = B A ± ■ B B ± where q± is perpendicular to the in- 
terface, x G [0, 1] is in fact an experimentally accessible 
quantity. This will be discussed below. 

The conductance in the normal phase is given by 

G N = I N /V = 2e 2 Kj2x p MxK T *P- Define r N(C) = 
Gn(C)/Gn(tt/2). Using Equation ©, we find r N (C) = 
1 + (1 - d) 2 /(l + d) 2 x cos C where d = N-/N+ is the ra- 
tio of the densities of states. The dependence on the an- 
gle C is similar to tunnelling magnetoresistance between 
ferromagnetsi^ and vanishes if A/+ = TV- . 

In the superconducting phase, the quasiparticle cur- 



rent for T 
where 



becomes I qp = ire Ea p -^a t a p 7a 



„ab 



7a A , B 



®(\eV\-(\ X $\ + \x*\)) eV^J 



1-5 



Xp 



(J) 
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K(m) and E(m) are the complete elliptic integrals of the 
first and second kind, respectively^ We have defined 
ox P = (Ixtl + \x B P \)/eV and 6 Xp = (|x£l " \x*\)/eV. 
This is a two-band generalisation of the one-band s-wave 
expression, 26 The usual one-band treshold at eV = x A + 
X B is replaced by at most four discontinuities. 

At zero voltage difference, the Josephson current be- 
comes 7j = 4ne £ Ap M A M B r Xp Y AB sin(tf B _ tfA) , w here 



1 Ap 



\xi\ + \xj 



K 



(lx£Mx B l) 2 
(Ixtl + lx B l) 2 



(8) 



£ = Jt/2 

- - c-n 




C-0 

£ = !l/2 



FIG. 1: (Color online) The quasiparticle current-voltage dia- 
gram at T = for angles C, = 0,n/2,n, x — 0.7 and d = 1. 
The graphs are displaced in the vertical direction for clarity. 
Left panel: |x— l/lx+l = 0.5. Right panel: \\- |/|x+ 1 — 0.01. 



This is a general two-band s-wave expression valid when 
interband Cooper pairs are absent. It is a straightforward 
generalisation of the standard one-band result dd 

We now consider the case of spin-orbit split bands. 
Usually, for equal systems one could let |x±| = |x±l = 
|x±|- This might not always be justified in our model, 
since the direction dependence of the gaps might depend 
on the nature of B^^) . However, let us again turn to the 

example above, where n A and n B point parallel to the in- 
terface. This makes k,p ~ q± equivalent directions even 
though n A ^ n B , at least in the isotropic approxima- 
tion. The above assumption should then be justified and 
will be used below. In addition, since the model (J2J con- 
tains interband pair scattering, we consider phase-locked 
bands where i9 + = + nn and n is zero or one. We do 
not investigate the possibility of small oscillations of the 
interband phase differenced 

Whenever |x+| = |x-l and A/+ = A/"_, / qp becomes 
independent of £ and equals the one-band result. In the 
case of unequal gaps, extra discontinuities should appear 
in the current-voltage characteristics. Figurcnshows the 
quasiparticle current at T = as function of voltage for 
the three angles £ = 0, tt/2 and it with x = 0.7 and d = 1. 
The graphs are displaced in the vertical direction for clar- 
ity. Note that as x — > all cases approach the £ = tt/2 
graph. In the left panel, |x-l/lx+l = 0.5. Discontinu- 
ities appear at eV = 2|x-|, (|x+l + |X-I) and 2|x+|. In 
the right panel, where |x-|/|x+l = 0.01, current can flow 
also for small voltages, although this almost vanishes as 
£ — > 7r. Smearing of the steps due to interband scatter- 
ing of quasiparticles should be negligible as long as the 
gap difference is well above fc^T. This breaks down in the 
limit of equal gaps, but then the current-voltage diagram 
collapse to the one-band result anyway. Anisotropic gaps, 
gap nodes and non-zero temperature may in general also 
lead to smearing. At non-zero temperatures, logarithmic 
singularities in 7 qp at eV = ||x+| — Ix-ll ma Y show up. 29 

Define |xm| = max(|x+|, |x-l) and |x m | = 
min (lx+l, |X-D- The ratio F = |Xm|/|XM| is easily found 
from the position of the first and last discontinuity. In 
addition, if Dx P denotes the jump in J qp at eV = (|xa| + 
\X P \), one finds d 2 F = D mm /D M M, where d = A^/A/'m- 
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These methods for finding F and d 2 F are independent 
of £ and x. Furthermore, x may be determined from 

-)/(Va; 



Da 



,D 



MM " 



D, 



xcosC = (2vA^Dmm 
for any angle 

The critical Josephson current 7j iC will 
also depend on (. There is a close anal- 
ogy to tunnelling magnetoresistance^ Define 
a = (-l) n 8j\f +J \f_T + _/iT(jSfl\x+\ + M^lX-l), where 
| a | = \a\(F,d) is monotonically increasing for F < 1, 
\a\(0,d) = and |a|(l,d) = 2d/(l + d 2 ) < 1. We find 
that a = — 1 results in Ij jC (tt/2) — 0. For a ^ — 1, we 
define 



rj(0 



_ MO 

" 'j,c (f ) 



1 



■ X COS ( 



(9) 



showing the possible modulation of the critical Josephson 
current with £. In addition, a = (1 — rj(0) + x)/(rj(0) — 
1 + x). The sign of a determines n and hence the relative 
sign between x+,q± an d X-,q±- One may then determine 
the ratio between the singlet and triplet components of the 
order parameters in a spin basis, since \rjq ± s\/\Vq± t| = 
(l + (-l)»F)/(l-(-l)«F). 

Note that if |x+| = |X-| and n = 0, rj(C) = r N (C). 
If in addition d = 1, Jj jC becomes independent of C and 
the one-band result^ is recovered. A modulation may 
be a result of unequal gaps (F ^ 1), unequal densities 
of states (d ^ 1) or both. In addition, no modulation 
of rj(£) could be interpreted both as x = and a = 1. 
Both these ambiguities should be distinguishable through 
the quasiparticle current-voltage characteristics. Consis- 



tency demands that \a\(F, d) found from the Josephson 
current fits F and d found from the quasiparticle current. 

To determine that the jumps in the quasiparticle cur- 
rent arise from spin-orbit split bands due to breakdown 
of inversion symmetry along a certain axis, several junc- 
tions with different relative orientations of those axes 
would be needed^ The synthesis and manipulation of 
such junctions thus represents a considerable experimen- 
tal challenge. However, building Josephson junctions 
with controllable crystallographic orientations was essen- 
tial to proving the cZ-wave symmetry of the order param- 
eter in the high-T c cuprates^ Also, the presence of two 
gaps may be possible to detect in other experiments, such 
as scanning tunnelling microscopy. 

In conclusion, we predict possible new effects in the 
tunnelling current between non-centrosymmetric super- 
conductors with significant spin-orbit splitting. Spin con- 
servation in the tunnelling barrier may then result in a 
modulation of the critical Josephson current when vary- 
ing the relative angle between the vectors describing ab- 
sence of inversion symmetry on each side. We have also 
shown that several discontinuities may appear in the 
quasiparticle current. We have argued that both these 
phenomena might help to determine the possibly exotic 
gap symmetry and pairing state of non-centrosymmetric 
superconductors. 
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